Cluster dynamical mean-field study of the Hubbard 
model on a 3D frustrated hyperkagome lattice 
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Abstract. We study the Hubbard model on a geometrically- frustrated hyperkagome lattice 
by a cluster extension of the dynamical mean field theory. We calculate the temperature (T) 
dependences of the specific heat (C) and the spin- lattice relaxation time (Ti) in correlated 
Ch 1 metallic region. C/T shows a peak at T = T p i and rapidly decreases as T — >• 0. On the other 

^ ' hand, 1/T±T has a peak at a higher temperature T P 2 than T p i, and largely decreases below T P 2, 

followed by the Korringa law 1/Ti oc T as T ^ 0. Both peak temperatures are suppressed and 
the peaks become sharper as electron correlation is increased. These behaviors originate from 
strong renormalization of the energy scales in the peculiar electronic structure in this frustrated 
I ■ system; a pseudo-gap like feature, the van-Hove singularity, and the flat band. The results are 

discussed in comparison with the experimental data in the hyperkagome material, Na^rsOs 
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1. Introduction 

How the electronic structure and thermodynamic properties evolve through the Mott transition 
has been a central issue in condensed matter physics [T]. Recently, insulating states close to the 
Mott transition have attracted considerable attention due to their exotic properties. Typical 
examples are found in a 2D organic conductor [2] and 3 He adsorbed on graphite [3j: Both are 
considered to be in the vicinity of the Mott transition having a tiny charge gap, and exhibit 
gapless spin-liquid behavior at low temperatures. If electron interaction is sufficiently strong 
and there is a large charge gap, a localized spin model gives a good description of the systems. 
Whereas, in the vicinity of the Mott transition, the expansion from the strong coupling limit 
suffers from slow convergence due to the small charge gap. In order to describe such systems, it 
rN \ is crucial to take account of charge fluctuations as well as spin ones. 

Among such insulators near the Mott transition showing exotic behaviors, we focus on 
Na4lr308 in this paper. The lattice structure of Na^^Os is similar to that of spinel oxides, 
while 1/4 of Ir 4+ pyrochlore sites are regularly replaced by non-magnetic Na + cations 0j. The 
Ir 4+ sublattice consists of a highly-frustrated three-dimensional (3D) network of corner-sharing 
triangles, which is called the 'hyperkagome' lattice [Fig. U[a)]. While the hyperkagome lattice 
shares a common local connectivity to the 2D kagome lattice, it has several distinct features 
such as the three dimensionality and structural chirality. It was reported that the charge gap 
of Na4lr308 is fairly small, suggesting that the compound is near the Mott transition: The 
magnitude of the gap is estimated at ~ 500 K 0]. 

Puzzling behaviors have been reported for the thermodynamic properties of Na^^Os 01- 
The Curie- Weiss temperature is large ~ 650 K, which is comparable to the magnitude of the 



charge gap. In spite of the large Curie- Weiss temperature, there is no clear sign of phase 
transition down to 2 K. C/T shows a peak at T ~ 20 K, and rapidly decreases as T —> 0. 
In contrast, the magnetic susceptibility takes a finite value as T —> 0, implying a macroscopic 
number of gapless magnetic excitations. As to the spin dynamics, the spin-lattice relaxation 
time, Ti, shows anomalous temperature dependence: 1/T\T shows a small peak at T ~ 250 
K, and decreases down to T ~ 50 K. Below T ~ 50 K, \jT\T shows a moderate temperature 
dependence, seemingly obeying the Korringa law [S]. 

In spite of several theoretical efforts to make sense of the wide breadth of experimental 
results, comprehensive understanding has not been reached yet. Instabilities toward spin nematic 
phase and Z2 spin liquid were discussed for the classical and quantum Heisenberg model on the 
hyperkagome lattice [HI HE]- Effects of the Dzyaloshinskii-Moriya interaction and anisotropic 
exchange interaction have also been studied [3]. 

A part of the difficulties faced by the theoretical studies stems from the fact that this 
compound locates in the vicinity of the Mott transition. The exchange interaction estimated from 
the Curie- Weiss temperature is comparable to the magnitude of the charge gap, which suggests 
an importance of including charge fluctuations in addition to spin ones. So far, however, all the 
theoretical analyses have been done for localized spin models. In this paper, in order to gain 
a new insight into this puzzling situation, we consider the problem from the opposite side, i.e., 
from the metallic region. We will explore how the electronic structure evolves as the electron 
correlation increases in the 3D frustrated hyperkagome systems. 

2. Model and Method 

We study the Hubbard model on the 3D frustrated hyperkagome lattice shown in Fig. [0(a) , 
whose Hamiltonian is given by 

H = ~t ^ (4a C jv + h - C + U 52 U it n il ~ V 52 W 
(i,j),(T i iv 

in the standard notations. We consider only the nearest-neighbor hopping and set t = 1 as 
an energy unit hereafter. The chemical potential \i is controlled so that the system is at the 
half filling (one electron per site on average). In Na^^Os, each Ir 4+ cation has one electron in 
three-fold t2 g orbitals on average. Our model may describe the case that the trigonal distortion 
splits the t2 g levels and the lower nondegenerate a\ g level plays a major role. 

To investigate the thermodynamics of this model, we adopt a cluster extension of the 
dynamical mean field theory jTU] . In this method, the original lattice problem is mapped to 
a cluster impurity problem. In the present case, in order to take account of the unique spatial 
geometry of the hyperkagome structure shown in Fig. Q](a), we consider a 12-sites cluster (the 
unit cell of the hyperkagome lattice), and apply the iterative perturbation theory (IPT) as an 
impurity solver. IPT is known to give a reliable description including the Mott transition for 
a single-impurity problem while it becomes poor for large U in the case of clusters. In 
the present study, we limit our investigations in the small to intermediate U regime where the 
method is more reliable. 

3. Results 

First let us consider the basic electronic structure of the hyperkagome model. In Fig. [U(b), we 
show the density of states (DOS) at U = in comparison with that of the 2D kagome case. 
The diverging peak at the upper band edge (£fl a t — 1.52) comes from four- fold-degenerate flat 
bands. In addition, two sharp peaks exist in the middle of the band, which are the van-Hove 
singularities at E v i ~ —0.48 and E V 2 — —2.48, respectively. These features are similar to the 
2D kagome case, however, there is an important difference: DOS of the hyperkagome model has 



a large dip near the Fermi energy (cj = 0), showing a pseudo-gap like feature. As a consequence, 
DOS at the Fermi energy is largely suppressed down to ~ 20% of that in the 2D kagome case. 
This small DOS might be relevant to consider the fact that Na^^Os is an insulator even though 
the electron correlation is expected to be weak for 5d electrons of Ir 4+ cations. 

When we switch on the electron correlation J7, these features of DOS near the Fermi energy 
undergo a strong renormalization by U. We show the results at low T in Fig. [2j Energies of 
the flat band and the van-Hove singularity, £fl at and E v i, are renormalized and approach the 
Fermi energy. These energy scales and their renormalization give rise to peculiar temperature 
dependences as we will see below. 





Figure 1. (a) Unit cell of the hyperkagome lattice. Figure 2. U dependence of the 
(b) Density of states of the noninteracting tight- density of states near the Fermi 
binding model on the 3D hyperkagome lattice in energy at T = 0.1. 
comparison with that on the 2D kagome lattice. 
Fermi energy is set at uj — 0. 



Next, we show the results for the specific heat (C) in Fig. [3l At U = 0, C/T has a peak 
at T = T p i ~ 0.3, and rapidly decreases for T < T p i. T p i correlates with the energy scale of 
the van-Hove singularity closer to the Fermi level, E Y \\ When T crosses |-E v i| 3 C shows a peak. 
The rapid decrease at lower T is ascribed to the small DOS at the Fermi level. For comparison, 
we show C/T for the 2D kagome model at U — in Fig. [3j which does not show such rapid 
decrease. As U increases, T p i shifts to lower T and the peak becomes sharper. As a result, the 
entropy is released more rapidly at low T for larger U. These interesting behaviors are closely 
related to the renormalization of DOS by U shown in Fig. [2j We note that the behaviors of C/T 
are similar to those observed in the experiment in Na^^Os |4J. 







Figure 3. Temperature dependence of Figure 4. Temperature dependence of 
C/T. The dashed line is the result at U = 1 /T\T. 
for the 2D kagome case for comparison. 



Finally, we discuss the spin relaxation time T\. We calculate \jT\T from the standard 
formula, 1/T\T oc lim^^o SqTr Imx H (q,^ + iO)/u)N. In order to obtain the transverse spin 
susceptibility, x H (q?^ + ^0), we approximate the irreducible vertex function in the Bethe- 
Salpeter equation by the RPA contribution. Figure [4] shows the results. 1/T\T shows a broad 
peak and decreases as T is lowered. At low T below ~ 0.05, it approaches a constant value, 
obeying the Korringa law. Although the peak shifts to lower T and the peak becomes sharper for 
larger U as seen in C/T, the peak temperature, T p 2, is different from T p i and higher than T p i. 
This might be attributed to the fact that \jT\T is more sensitive to the flat band at a higher 
energy whose singularity is much stronger than the van-Hove singularity. Note that 1/T\T is 
proportional to the square of DOS at U = 0, in contrast to C/T which is proportional to DOS 
itself. Our results show similar features to \/T\T of Na^^Os such as a peak at high T followed 
by a large drop and a Korringa-like behavior at the lowest T [5j. 

4. Discussion and Concluding Remarks 

We have obtained characteristic temperature dependences of C/T and \ jT\T\ Both show a peak 
followed by a rapid decrease at low T, and the peaks shift to lower T as increasing U. These 
behaviors are explained by the strong renormalization effect by U on the peculiar structure of 
DOS in the hyperkagome model. Particularly important factors are the pseudo-gap like feature 
and two singularities, the van-Hove singularity and the flat band. 

It is rather surprising that these results in the metallic regime are qualitatively similar to 
those observed in the hyperkagome material, Na^^Os, which is an insulator with a small charge 
gap. The agreement implies a possibility that the renormalized singularities of DOS survive and 
affect the thermodynamic properties even in the insulating regime when the system is in the 
vicinity of the Mott transition. Note that although the Mott transition itself is known to be 
of first order, it is like a liquid-gas transition where the finite-temperature properties can be 
smoothly connected by bypassing the critical endpoint [I]. 

To explore this fascinating scenario, it is necessary to extend our analysis to the insulating 
region beyond the Mott transition. Such study by employing an improved impurity solver is 
in progress. There, it is interesting to examine whether the vanishing or very tiny C/T, the 
Korringa-like behavior of 1/TiT, and a finite value of the magnetic susceptibility at T are 
comprehensively explained. Considering the real compound, it is also interesting to consider 
effects of orbital degrees of freedom, spin-orbit interaction, and disorder due to Na deficiency. 
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